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Abstract:  A  theoretical  study  has  been  made  of  the  shallow  donor  and  Wannier 
exciton  within  a  one-dimensional  quantum  well.  The  variational  method  was 
used  with  a  cylindrical  oaussian  basis,  sc..  It.  order  to  facilitate  comparison 
with  future  experimental  measurements  of  excited  states  of  these  systems,  an 
external  magnetic  field  was  assuaed  perpendicular  to  the  interfaces  between 
the  barrier  material  and  the  well.  Calculations  reveal  that  the  choice  of 
matching  conditions  used  at  the  interfaces  has  little  effect  on  the  binding 
energies  of  the  ground  or  first  few  excited  states  of  the  shallow  donor,  ex¬ 
cept  for  well  widths  considerably  smaller  than  the  effective  Bohr  radius. 

The  results  of  calculations  of  the  shallow  donor  ground  and  first  few  excited 
states  are  presented  for  a  variety  of  well  sizes  and  magnetic  field  strengths. 
Similar  results  are  given  for  the  ground  state  of  the  Wennler  exciton. 
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Final  Report 


Abstract:  A  theoretical  study  has  been  made  of  the  shallow  donor  and  Wannier 
exciton  within  a  one-dimensional  quantum  well.  The  variational  method  was 
used  with  a  cylindrical  Gaussian  basis  set.  In  order  to  facilitate  comparison 
with  future  experimental  measurements  of  excited  states  of  these  systems,  an 
external  magnetic  field  was  assumed  perpendicular  to  the  interfaces  between 
the  barrier  material  and  the  well.  Calculations  reveal  that  the  choice  of 
matching  conditions  used  at  the  interfaces  has  little  effect  on  the  binding 
energies  of  the  ground  or  first  few  excited  states  of  the  shallow  donor,  ex¬ 
cept  for  well  widths  considerably  smaller  than  the  effective  Bohr  radius. 

The  results  of  calculations  of  the  shallow  donor  ground  and  first  few  excited 
states  are  presented  for  a  variety  of  well  sizes  and  magnetic  field  strengths. 
Similar  results  are  given  for  the  ground  state  of  the  Wannier  exciton. 
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I. 


Introduction 


Recent  advances  in  epitaxial  crystal  growth  techniques  have  made  it  poss¬ 
ible  to  grow  systems  consisting  of  alternate  layers  of  two  different  semicon¬ 
ductors,  with  sharp  interfaces.  One  of  the  most  common  of  these  layered  systems 
consists  of  GaAs  alternated  with  Ga,  Al  As.  The  difference  in  band  gap  of  the 

i.  "“X  X 

two  semiconductors  leads  to  conduction  and  valence  band  discontinuities  at  the 
interfaces.  This  effectively  produces  a  periodic  square-well  potential  with 
period  controllable  during  growth.  This  large  scale  periodic  potential  results 
in  the  breaking  the  usual  electron  energy  bands  into  subbands,  and  has  signi¬ 
ficant  effects  on  the  transport  properties  and  bound,  atomic-like  systems  of 
these  materials.  The  research  reported  herein  deals  with  the  effects  of  quan¬ 
tum  wells  upon  hydrogen-like  systems,  specifically  shallow  donors  and  Wannier 
excitons. 

The  study  of  such  systems  is  still  in  its  infancy,  rhe  first  calculations 
of  the  ground  state  binding  energy  of  shallow  impurities  in  quantum  wells  as- 

1  2 

sumed  an  infinite  barrier  between  the  well  material  and  the  barrier  material.  * 

3  4  5 

Later  calculations  by  Mailhiot  e£  al.  *  and  by  Greene  and  Bajaj  replaced  the 
infinite  barriers  by  finite  band  discontinuities  dependent  upon  x,  the  fraction¬ 
al  Al  content,  and  determined  energies  for  the  first  few  excited  states  as  well 
as  the  ground  state.  The  latter  authors  in  particular  found  significant  quali¬ 
tative  and  quantitative  differences  between  the  behavior  of  2s-  and  2p+-like 
excited  states  compared  to  the  2po-like  state. 

Assuming  infinite  quantum  wells.  Bastard  £t  £l.^  performed  similar  calcul¬ 
ations  of  the  binding  energy  of  the  ground  state  of  Wannier  excitons  in  quantum 
well  structures.  Shortly  afterward,  Greene  and  Bajaj^  reported  results  of  the 
exciton  ground  state  in  finite-well  GaAs-Ga^^Al^As  structures,  and  very  re- 

O 

cently  extended  the  work  to  include  a  few  low-lying  excited  states. 
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One  of  the  major  reasons  for  calculating  the  binding  energies  of  these 
hydrogen-like  systems  is  the  expectation  that  comparison  between  theory  and 
experiment  can  give  more  information  about  the  general  quality  of  the  layered 
materials  and  about  the”  nature  of  the  interface  between  the  two  semiconductors. 
However,  it  is  difficult  to  experimentally  determine  donor  and  exciton  excited 
state  energies  due  to  the  large  line  widths  of  ground-to-excited-state  trans¬ 
itions.  In  bulk  GaAs  these  line  widths  can  be  reduced  by  the  application  of 

q 

an  external  magnetic  field.  It  was  the  anticipation  that  a  similar  technique 
would  be  used  in  these  quantum  well  structures  that  led  to  the  investigation  of 
the  Zeeman  levels  of  shallow  donors  and  Wannier  exciton  in  quantum  wells. 

The  objectives  of  the  research  were  the  following: 

1)  To  study  how  the  treatment  of  the  different  semiconductor 
effective  masses  and  dielectric  constants  affects  the  binding 
energies  of  shallow  donor  states, 

2)  To  investigate  the  possibility  of  using  the  energy  difference 
between  the  shallow  donor  2p+-like  and  2po~like  states  to  pro¬ 
vide  an  independent  determination  of  the  conduction  band  dis¬ 
continuity, 

3)  To  determine  the  binding  energies  of  the  ground  and  several 
excited  states  of  the  shallow  donor  at  the  center  of  the 
GaAs-Ga^  xAlxAs  quantum  well,  as  a  function  of  well  size  and 
applied  magnetic  field, 

A)  To  obtain  results  for  the  Wannier  exciton  ground  state  energy 
in  the  GaAs-Gaj,  xAlxAs  quantum  well  as  a  function  of  well 
size  and  applied  magnetic  field. 

The  results  are  given  in  Section  III,  and  specific  accomplishments  are 
summarized  in  Section  IV  of  this  report. 
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II.  Theory  and  Method  Used 


A.  Shallow  Donor 

The  effective  mass  approximation  is  assumed  valid  for  the  GaAs-Ga^^Al^As 

quantum  wells  under  consideration.  With  this  approximation  the  Hamiltonian 

for  a  hydrogenic  donor  in  a  magnetic  field,  applied  in  a  direction  prependicular 

to  the  layers,  can  be  written 

H*— — r  (-itiV  -  —  )2  -  ~~  +  V'w(z)  ,  (1) 

2m*  c  eQr  * 

The  donor  impurity  is  taken  to  be  located  at  the  center  of  the  well.  The  well 
is  a  finite  square  well 


^0 


|z  |  <  L/2 

|z  |  >  L/2  , 


(2) 


where  L  is  the  width  of  the  quantum  well,  and  VQ  is  the  conduction  band  discon¬ 
tinuity.  Using  the  cylindrical  gauge,  where 

l  4  -> 

A*  y  B  x  r. 


and  a  cylindrical  coordinate  system  with  Z-axis  along  the  magnetic  field  direc¬ 
tion,  yields 

H-  -  ^V2  -  |  +  YLz  +  i  YV  +  Vw(z)  ,  (3) 

where 

r-  (p2  +  z2)2 

and  L  is  the  z-component  of  the  angular  momentum  operator  (in  units  of  ti) . 

X 

The  Hamiltonian  in  Gq.  (3)  is  written  in  units  of  the  effective  Rydberg  in  bulk 


GaAs,  m*e% 

R*- 

and  the  unit  of  length  is  the  GaAs  effective  Bohr  radius, 

Veo 

a*-  m"*e5 

The  mfi  in  Eq.  (3)  is  the  effective  mass  of  the  electron  in  units  of  the  bulk 

GaAs  electron  effective  mass,  m*.  For  the  GaAs-Ga.  A1  As  structures,  m  >1 

1-x  x  e 

within  the  well,and  mg>l  (depending  on  x)  in  the  barrier  material.  Finally, 
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the  parameter y is  a  dimensionless  measure  of  the  magnetic  field,  given  by 

Y.  ehB _ 

1  2m*cR* 

The  conduction  band  discontinuity  (or  barrier  height) in  Eq.  (2)  is  assumed 

to  be  85%  of  the  total  band  gap  discontinuity  determined  from  the  empirical 
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expression  AE  =  1# 155x  +  .37X  eV.  (4) 

O 

For  typical  values  of  x  and  L,  the  barrier  height  is  considerably  larger  than 
the  effective  Rydberg.  Thus,  the  energy  associated  with  the  Coulomb  inter¬ 
action  will,  except  for  large  well  widths,  be  small  compared  to  the  subband 
energy  due  to  the  square  well.  Because  of  this,  it  is  helpful  to  explicitly 
factor  the  solution  to  the  ground  state  of  an  electron  in  the  one-dimensional 
square  well  out  of  the  donor  wave  function  V: 

V-  f(z)G(p,*,<fr)  .  (5) 


The  (unnormalized)  square  well  solution  is  given  by 
Ceos  kz  |z  I  <  L/  2 

£(z)“  (Ae"K  lz  •  |Z1  *  L/2  (6) 

The  parameter  k  is  determined  from  the  energy  of  the  first  subband,  and  the 

constants  A  and  <  are  fixed  by  the  matching  conditions  at  the  interface.  It 

is  assumed  that  F(z)  and  (1/m  )  3f/3z  are  continuous  across  the  interface.  As 

discussed  in  Sectionlll,  however,  the  results  are  not  very  sensitive  to  the 

matching  conditions  chosen. 

For  wells  less  than  several  effective  Bohr  radii  in  width,  the  factoriza¬ 
tion  in  Eq. (5)  is  particularly  appropriate  because  it  enables  the  removal  of 
the  relatively  large  square  well  energy  that  otherwise  would  numerically  over¬ 
whelm  the  smaller  Coulomb  energy  contribution.  For  large  well  widths,  where 
the  energy  of  the  subband  is  less  than  or  of  the  same  order  as  the  Coulomb 
contribution,  the  factorization  is  of  little  or  no  value  because  higher  sub¬ 
bands  can  be  mixed  with  the  lowest  subband  by  the  Coulomb  potential.  However 
this  does  not  cause  any  significant  problem  provided  the  function  G(p,z,<f>) 
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has  sufficient  variational  flexibility. 

Since  the  Hamiltonian  of  £q.(3)  is  cylindr ically  symmetric,  the  Z-compon- 
ent  of  the  angular  momentum  is  a  good  quantum  number.  The  dependence  of  the 
wave  function  thus  has  the  form  exp(im<i),  where  m  is  an  integer.  If  the  donor 
impurity  is  located  at  the  center  of  the  well,  the  Hamiltonian  is  also  invariant 
under  reflection  through  the  origin.  The  wave  function  has  a  definite  parity. 
Using  this  knowledge,  the  function  G(p,z,$)  can  be  written  in  the  form 
G(p,z,<J>)“  p!mleim4,zq  l  A^  G11(p,z)  .  (7) 

i.j 

The  parity  of  this  function  is  determined  by  m+q ,  where  q  »  0,  1. 

The  basis  functions  (p,  z)  are  taken  to  be  products  of  Gausslans  in 
P  and  z  variables:  G^fp.z)-  e'(6+  <H>P2  .  (8) 

This  choice  was  made  because  of  the  success  of  a  similar  basis  set  which  Aldrich 

and  Greene H  applied  to  the  problem  of  hydrogen  in  a  uniform  magnetic  field. 

They  found  that  the  Gaussian  basis  set  yielded  good  results  for  the  ground  and 

several  excited  states  throughout  the  range  0  £  y  <.  10  . 

The  set  of  parameters }  used  in  this  work  is  given  in  Table  1.  They  are 

12 

taken  from  the  results  of  Huzinaga,  who  did  a  detailed  study  of  the  use  of 
Gaussian  basis  functions  in  the  calculation  of  electron  energies  in  hydrogen. 
This  set  of  {  }  gives  energies  for  the  Is,  2s  and  2p  f ree-hydrogen  states 

accurate  to  within  0.001  Rydbergs.  The  parameter  3  was  varied  in  each  case  to 
minimize  the  energy.  It  is  primarily  determined  by  the  size  of  the  magnetic 
field.  Calculations  were  also  made  using  3  as  a  multiplicative  variational 
parameter  in  the  argument  of  the  p  Gaussian,  rather  than  the  additive  one  of 
Eq.(8).  The  results  were  poorer  (gave  smaller  binding  energies)  by  up  to  ten 
per  cent,  particularly  for  larger  magnetic  field  strengths.  The  additive  $ 
differentially  distorts  the  individual  p  Caussians,  in  a  way  consistent  with 
the  physics  of  the  magnetic  field  (which  affects  the  bound  electron  least 
strongly  when  it  is  near  the  donor).  The  multiplicative  3  ,  by  contrast. 


provides  the  same  scaling  for  each  of  the  p  Gaussians. 


The  results  presented  in  the  next  section  used  A^j*0  for  i^j ;  the  energie 
were  then  determined  by  solving  the  matrix  eigenvalue  equation. 


Hl|>  -  Elty  , 

where  H  and  U  are  the  Hamiltonian  and  overlap  matrices.  For  m*q*0  (even  parity), 
the  number  of  basis  functions  (and  hence  the  order  of  the  Hamiltonian  and  over¬ 
lap  matrices)  is  five.  For  m*l,  q=0  and  m=0,  q=l  (odd  parity),  there  are  four 
basis  functions. 

A  number  of  runs  at  various  values  of  y  and  L  were  also  made  with  A^j=0 
for  ji— j ( > 1  and  non-zero  otherwise.  The  energies  were  obtained  as  before.  In 
this  case  there  are  13  and  10  basis  functions  for  the  even  and  odd  parity  cases, 
respectively.  The  larger  number  of  basis  functions  yielded  improvements  in  the 
calculated  binding  energies  of  about  0.001  effective  Rydbergs  or  less,  which 
gives  confidence  in  the  accuracy  of  the  results. 

B.  Wannier  Exciton 

If  the  quantum  well  associated  with  the  valence  band  discontinuity  is 
assumed  large  enough  to  split  the  degeneracy  of  the  GaAs  valence  band,  two  exci- 
tons  (known  as  the  light-hole  and  heavy-hole  excitons)  arise.  In  each  case  the 
reduction  of  the  Hamiltonian  follows  pretty  much  the  same  steps  as  for  the  donor. 
The  major  difference  arises  because  the  hole  moves  in  the  Coulomb  potential  and 
its  own  quantum  well  with  an  effective  mass  which  is  not  isotropic.  The  quan¬ 
tity  is  taken  to  be  the  heavy  (+)  or  light  (-)  hole  effective  mass  in  the 


direction  perpendicular  to  the  layers  (Z-axis) ,  and  is  the  electron-hole  re¬ 
duced  effective  mass  corresponding  to  the  heavy  (+)  or  light  (-)  hole  bands  in 

«-he  plane  perpendicular  to  the  Z-axis.  These  masses  are  given  in  terms  of  the 
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Kohn-Luttinger  band  parameters  y  ^  and  as 


and 


—  -  —  +  <Y,  ±  Y  ) 
"Vi  12 


(9a) 


l 

m 


Y,  +  2y« 


The  masses  in  these  equations  are  in  units  of  the  free  electron  mass. 


In  terms  of  these  masses,  the  Hamiltonian  for  the  exciton  in  the  quantum 


wells  formed  by  the  conduction  and  valence  band  discontinuities  is  given  by 

2 

±  92 

(10) 


„  1  3  3  .  1  32  “I  u±  32  u± 

h’"Lp^pp3"p  P2  W  J  '  ^7z l  ~  ~ 
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-  -  +  yL  +  fY2P2  +  V  (z  )  +  V.  (z,  )  . 

r  1  z  4  ew  e  hw  *h 


The  symbols  p,  <j> ,  and  r  are  electron-hole  relative  coordinates,  and  the 
Hamiltonian  is  expressed  in  units  of  the  effective  Rydberg  defined  with  elec¬ 
tron-hole  reduced  mass,  u+. 

Since  the  exciton  Hamiltonian  of  Eq.  (10)  has  cylindrical  symmetry  (as 
does  that  for  the  donor,  Eq.  (3))  the  same  type  of  variational  solution  was 
chosen.  The  only  difference  is  that  the  ground  state  solution  for  the  hole  in 
a  one-dimensional  square  well  was  factored  out  of  the  variational  trial  funct¬ 
ion,  in  addition  to  that  for  the  electron: 


f^re’rh^  =  fe^Ze^fh^Zh^  ^P*2*^* 


(ID 


The  function  G(p,  z,  d>>  >  which  is  assumed  to  be  a  function  only  of  electron- 
hole  relative  coordinates,  is  of  the  same  form  as  Eq. (7)  and  (8).  For  the  re¬ 
sults  of  the  next  section  A,  ,«0  for  i^i,  and  the  are  determined  from  the 

i  j  J  ii 

solution  of  the  eigenvalue  problem. 

The  set  of  parameters  {ct^}  Is  the  same  as  before  since  the  calculation 
is  done  using  effective  atomic  units.  However,  the  Hamiltonian  and  overlap 
matrices  now  involve  integrals  over  both  electron  and  hole  positions. 

III.  Results  and  discussion 

To  examine  the  importance  of  the  different  effective  masses  of  GaAs  and 
Ga^  xAlxAs,  two  separate  sets  of  calculations  of  the  donor  binding  energy  of  the 
ground  and  several  excited  states  were  performed  for  x*.15,  .30;  Y  *0.0, 


a'jiS 


-\  \  A 

i  *_a  \j* 
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.1,  .2,  .3,  .5,  .75,  1.0,  2.0,  5.0;  and  various  values  of  L  from  50-1000A. 

One  set  of  calculations  used  the  GaAs  value  for  the  electron  effective  mass  in 
both  the  barrier  and  well  semiconductors.  The  interface  matching  conditions 
then  required  continuity  of  the  function  f  of  Eq.  (6)  and  3f/3z  across  the  inter 
face.  These  are  the  conditions  assumed  by  Greene  and  Bajaj.^  The  other  set  of 
calculations  used  the  following  expression  for  the  Ga^  ^Al^As  effective  mass 

4 

(m  ,  in  units  of  the  free  electron  mass) : 
e 

m  -  0.067  +  0.83x  (11) 

e 

The  matching  conditions  at  the  interface  in  this  case  were  that  f(z)  and  (1/m  ) 

e 

3f/3z  are  continuous  across  the  interface.  These  conditions  are  similar  to 

3  4 

those  used  by  Mailhiot  et  al.  *  and  have  been  justified  for  the  GaAs-Ga^^Al^As 

14 

system  by  the  band  calculations  of  Ando  and  Mori. 

The  differences  between  the  binding  energies  obtained  in  the  two  sets  of 
calculations  are  very  small.  In  most  cases  they  differ  by  less  than  1%.  The 
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worst  case  occurs  for  x*.15,  L*50A  with  small  Y  .  The  differences  there  are 
about  3%  for  the  ls-like  and  2 p+-like  states.  Thus  it  appears  that  as  far  as 
the  binding  energy  of  the  shallow  donor  is  concerned,  it  does  not  matter  much 
what  matching  conditions  are  used,  except  for  very  small  L. 

The  original  proposal  for  this  work  indicated  that  a  study  of  the  effects 
of  the  different  dielectric  constants  of  the  two  semiconductors  would  be  under¬ 
taken.  This  intention  was  prompted  by  the  incorrect  treatment  of  the  potential 

4 

with  the  two  dielectric  constants  in  the  early  work  of  Mailhiot  £t  al.  After 
the  proposal  was  submitted,  Maihiot  et.  al .  published  the  results  of  more  com¬ 
plete  calculations^  which  corrected  the  error.  Because  their  results  are  very 
close  to  the  results  of  Greene  and  Bajaj^  who  used  the  same  dielectric  constant 
for  both  GaAs  and  Ga^  ^Al^As,  it  was  decided  that  additional  calculations  using 
different  dielectric  constants  would  not  be  worthwhile. 


’  v  '  O*  1.  i  . 
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With  regard  to  the  second  research  objective,  a  number  of  calculations  of 
the  binding  energy  of  the  2p+-like  and  2pQ-like  states  were  made,  varying  the 
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barrier  height.  Two  different  well  widths,  L=500A  and  L=750A  were  used.  The 

conduction  band  discontinuity  was  varied  in  5%  increments  from  70%  to  90%  of 

the  total  band  gap  discontinuity,  AE  .  The  hope  was  that  because  of  the  signi- 
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ficantly  different  behavior  of  these  two  states  as  a  function  of  L,  they  could 

be  used  in  conjunction  with  experiment  to  accurately  determine  the  conduction 

band  discontinuity.  The  results  of  these  calculations  indicate  that  this  is 

likely  not  the  case.  The  energy  difference  between  the  two  states  varies  rather 

slowly  with  barrier  height;  the  change  form  .7AE  to  .9AE  is  only  0.01  R*,  or 

g  g 

0.058  meV.  This  small  variation  is  not  enough  larger  than  the  expected  ac¬ 
curacy  of  the  calculations  to  be  acceptable. 

Tables  2  -  7  give  the  results  of  calculations  of  the  binding  energies  of 

the  ls-like,  2p+-like,  and  2pQ-like  shallow  donor  states  for  a  variety  of  well 
widths  (L)  and  magnetic  field  strengths  (y).  Two  values  of  A1  fractions,  x*.15 
and  x»«.30,  are  tabulated.  Tables  8-11  give  similar  results  for  the  ground 
state  binding  energies  of  the  light-and  heavy-hole  exciton.  Accuracy  of  all 
results  is  estimated  to  be  about  0.01  effective  Rydbergs. 


< 


IV.  Sunmary 

The  major  accomplishments  of  the  research  are  listed  below: 

1.  The  binding  energy  of  the  donor  ground  and  first  few 
excited  states  has  been  shown  to  be  very  insensitive 
to  the  interface  matching  conditions  used  to  take  in¬ 
to  account  the  affective  mass  difference  between 

electrons  in  GaAs  and  those  in  Ga,  A1  As,  at  least 

1— x  x 

for  well  widths  between  50  and  1000  A. 

2.  The  energy  difference  between  the  2p±-like  state  and 
the  2pQ-like  state  is  too  weakly  dependent  on  the  con¬ 
duction  band  discontinuity  to  use  that  difference  as 

a  diagnostic. 

3.  The  binding  energies  of  the  shallow  donor  ground  and 
first  few  excited  states  have  been  tabulated  for  a 
variety  of  well  widths  and  heights  and  magnetic  field 
strengths. 

4.  The  ground  state  binding  energies  of  a  Wannier  exciton 
in  a  quantum  well  has  similarly  been  tabulated  for 
various  well  widths  and  heights  and  magnetic  field 
strengths. 

The  work  was  performed  by  the  principal  investigator.  Dr.  Ronald  L.  Greene 
in  collaboration  with  Dr.  K.  K.  Bajaj  of  the  Air  Force  Avionics  Laboratory.  Two 
manuscripts  reporting  the  results  of  these  investigations  are  currently  be 
written  and  will  be  submitted  for  publication  in  Physical  Review  B. 
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Units  are  GaAs  effective  Rydbergs  (5.81  meV).  The  parameter  y  is  proportional  to  the  magnetic  field 


